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Electron, muon and tau magnetic moments: a theoretical update 
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a Dipartimento di Fisica, Universita di Padova and INFN, Via Marzolo 8, 35131 Padova, Italy 

Recent Standard Model predictions for the anomalous magnetic moments of the electron, muon and r lepton 
are reviewed and compared to the latest experimental values. 



1. INTRODUCTION 

The gyromagnetic factor g is defined by the re- 
lation between the particle's spin s and its mag- 
netic moment /2, 



e _ 
2m 



(1) 



where e and m are the charge and mass of the 
particle. In the Dirac theory of a charged point- 
like spin-1/2 particle, g = 2. Quantum Electro- 
dynamics (QED) predicts deviations from Dirac's 
value, as the charged particle can emit and reab- 
sorb virtual photons. These QED effects slightly 
increase the g value. It is conventional to express 
the difference of g from 2 in terms of the value 
of the so-called anomalous magnetic moment, a 
dimensionless quantity defined as a = (g — 2)/2. 

The measurements of the anomalous magnetic 
moment of the electron and the muon recently 
reached the fabulous relative precision of 0.7 parts 
per billion (ppb) [1] and 0.5 parts per million 
(ppm) [2,3], respectively. The theoretical predic- 
tion of the electron one, a e , is rather insensitive 
to strong and weak interactions, hence providing 
a stringent test of QED and leading to the most 
precise determination of the fine-structure con- 
stant a to date [4,5]. On the other hand, the 
muon anomalous magnetic moment allows to 
test the entire Standard Model (SM) and scruti- 
nize viable alternatives to this theory, as each of 
its sectors contributes in a significant way to the 
total prediction. Compared with a e , is also 
much better suited to unveil or constrain "New 
Physics" (NP) effects [6]. Indeed, for a lepton 
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/, their contribution to a; is generally expected 
to be proportional to mf/A , where m; is the 
mass of the lepton and A is the scale of NP, thus 
leading to an (m^/m e ) 2 ~ 4 x 10 4 relative en- 
hancement of the sensitivity of the muon versus 
the electron anomalous magnetic moment. This 
more than compensates the much higher accuracy 
with which the g factor of the latter is known. 
The anomalous magnetic moment of the r lep- 
ton, a T , would suit even better; however, its di- 
rect experimental measurement is prevented by 
the relatively short lifetime of this lepton, at least 
at present. The existing limits are based on the 
measurements of the total and differential cross 



sections of the reactions e + e 



and 



e + e — > Z — > t + t 7 at LEP energies. The most 
stringent limit, -0.052 < a^ xp < 0.013 at 95% 
confidence level, was set by the DELPHI collabo- 
ration [7] , and is still more than an order of mag- 
nitude worse than that required to determine a T . 

The SM prediction af M (i = e, ju or r) is usually 
split into three parts: QED, electroweak (EW) 
and hadronic. Here we provide a summary of the 
present status of these contributions, and a com- 
parison with the current experimental values. 

2. QED EFFECTS: PREAMBLE 

The QED part of the anomalous magnetic mo- 
ment ai of a charged lepton I = e, /i or t arises 
from the subset of SM diagrams containing only 
leptons and photons. For each of the three lep- 
tons /, of mass mi, this dimensionless quantity 
can be cast in the general form [8] 



QED_ 



mi 



A 1 +A 2 [ — ) + A 2 [ — 



mi 
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where mj and m k are the masses of the other 
two leptons. The term Ai, arising from diagrams 
containing photons and leptons of only one fla- 
vor, is mass and flavor independent. In contrast, 
the terms A 2 and A 3 are functions of the indi- 
cated mass ratios, and are generated by graphs 
containing also leptons of flavors different from 
I. The contribution of a lepton j to af EO is sup- 
pressed by (mf/rrij) if m 3 »m;, while it contains 
a logarithmic enhancement factor \n(mi/nij) if 
rrij <C mi. The muon contribution to a^ ED , for 
example, is thus power suppressed by a factor 
(ml/m 2 ^) ~ 10~ 5 ; nonetheless, this effect is much 
larger than the tiny experimental uncertainty on 
a e [1] (see Sec. 3). On the contrary, the QED 
parts of a^.T beyond one-loop are dominated by 
the mass-dependent terms. 

The functions Ai (i = l,2,3) can be expanded 
as power series in u/tt and computed order- by- 
order 

*_^(H) + ^(2)\^(2)V... (3) 

Only one diagram is involved in the evaluation 
of the lowest-order (first-order in a, second-order 
in the electric charge) contribution — it provides 

(2) 

the famous result by Schwinger A\ ' = 1/2 [9j. 
The mass-dependent coefficients A 2 and A 3 are of 
higher order. All results discussed below were de- 
rived using the latest CODATA [10] recommended 
mass ratios: m e /m^ = 4.836 33167(13) x 1(T 3 , 
m e /m T = 2.875 64(47) x 1(T 4 , m^/m, = 
206.768 2838(54), m p /m T = 5.945 92(97) x 10~ 2 , 
m T /m e = 3477.48(57), m T /m^ = 16.8183(27). 
The value for m T adopted by CODATA in Rcf. [10] 
(m T = 1776.99 (29) McV) is based on the 2002 
PDG result [11]. This PDG result remains un- 
changed to date [12,13]. 

3. ELECTRON 

3.1. QED Contributions 

Seven diagrams contribute to the fourth-order 
coefficient , one to A^ {m e /m^) and one to 
A 2 (m e /m T ). As there are no two-loop diagrams 
contributing to a® ED that contain both virtual 
muons and taus, ^'(rrie/m^, m e /m T ) = 0. The 
mass-independent coefficient has been known for 



almost fifty years [14]: 

A[ 4) = -0.328 478 965 579 193 78 ... . (4) 

The coefficient of the two-loop mass-dependent 
contribution to af ED , A^\l/x), with x = mj/mi, 
is generated by the diagram with a vacuum po- 
larization subgraph containing the virtual lepton 
j. This coefficient was first computed in the late 
1950s for the muon g— 2 with x = m e /m fl ^_l 7 ne- 
glecting terms of 0(x) [15]. The exact expression 
for < x < 1 was reported by Elcnd in 1966 [16]. 
However, for many years its numerical evaluation 
was considered tricky because of large cancella- 
tions and difficulties in the estimate of the ac- 
curacy of the results, so that common practice 
was to use series expansions instead [17-19]. Tak- 
ing advantage of the properties of the dilogarithm 
Li 2 (z) = -f*(dt/t) ln(l - t) [20], the exact result 
was cast in [21] in a very simple and compact 
analytic form, valid, contrary to the one in [16], 
also for x > 1 (the case relevant to aJ? ED and part 
of a^ ED ). Its numerical evaluation with the mass 
ratios given in Sec. 2 yields [5] 

AfK/mJ = 5.197 386 70(28) x 10~ 7 (5) 

A {4) {m e /m T ) = 1.837 62(60) x 10" 9 , (6) 

where the errors are only due to the uncertainties 
of the mass ratios. The results of Eqs. (5) and (6) 
are equal to those obtained with a series expan- 
sion in powers of the mass ratio y and lny, with 

y«i [io]. 

Adding up Eqs. (4), (5) and (6) one gets the 
two-loop QED coefficient [5] 

Cf) - 4 4) +4 4) (m e /TO M )+4 4) (m e /m T ) 
= -0.328 478 444 002 90(60). (7) 

The mass-dependent part of is small but not 
negligible, giving a relative contribution to the 
theoretical prediction of the electron g— 2 of 2.4 
ppb. This value is much larger than the 0.7 ppb 
relative uncertainty very recently achieved in the 
measurement of a e [1]. The uncertainties in 

A^ {m e /m^) and A^ (m e /m T ) are dominated by 
those in the latter and were added in quadrature. 
The resulting error Sci 4) = 6xl0~ 13 leads to a to- 
tally negligible O(10~ 18 ) uncertainty in the a^ ED 
prediction. 
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More than one hundred diagrams are involved 
in the evaluation of the three-loop (sixth-order) 
QED contribution. Their analytic computation 
required approximately three decades, ending in 
the late 1990s. The coefficient A^' arises from 72 
diagrams. Its exact expression is mainly due to 
Remiddi and his collaborators [22,23]; its numer- 
ical value is 

Af ] = 1.181241456 587..., (8) 

in very good agreement with previous results ob- 
tained with numerical methods [24] . 

The calculation of the exact expression for the 
coefficient (mi /rnj) for arbitrary values of the 
mass ratio mi/rrij was completed in 1993 by La- 
porta and Remiddi [25,26] (earlier works include 
Refs. [27]). Let us focus on a« ED (l = e, j = H,r). 
This coefficient can be further split into two parts: 
the first one, A^ (mi jvcij , vac) , receives contribu- 
tions from 36 diagrams containing either muon 
or t vacuum polarization loops [25], whereas the 
second one, A^ (mi /mj , lbl) , is due to 12 light- 
by-light scattering diagrams with either muon or 
t loops [26]. The exact expressions for these co- 
efficients are rather complicated, containing hun- 
dreds of polylogarithmic functions up to fifth de- 
gree (for the light-by-light diagrams) and complex 
arguments (for the vacuum polarization ones) — 
they also involve harmonic polylogarithms [28]. 
Series expansions were provided in Ref. [26] for 
the cases of physical relevance. 

Using the recommended mass ratios given in 
Sec. 2, the following values were recently com- 
puted from the full analytic expressions [5] : 

4 6) (m e /m M ) = -7.373 94164(29) x 10" 6 (9) 
A^\m e /m T ) = -6.5819(19) x 10" 8 . (10) 

Equations (9) and (10) provide the first evalua- 
tion of the full analytic expressions for these coef- 
ficients with the CODATA mass ratios of [10]; they 
are almost identical to the results A\ '(m e /m^) = 
-7.373 94158(28) x 10~ 6 and A^ ] (m e /m T ) = 
-6.5819(19) x 10~ 8 obtained in [10] via the ap- 
proximate series expansions in the mass ratios. 
The small difference between A^\m e /m^) of [10] 
and Eq. (9) mainly origins from the 0((m e /m p ) 6 ) 



term in the series expansion of A ( 2 (m e /m^, lbl); 
indeed, due to its smallncss, this term was ne- 
glected in the expansions of Ref. [26] used in [10]. 
The 0((m e /m,j,) 6 ) and 0((m e /m M ) 8 ) terms were 
explicitly provided in Ref. [5] expanding the ex- 
act Laporta-Rcmiddi expression for the sum of 
light-by-light and vacuum polarization contribu- 
tions for mi/mj <C 1 (see also [25,29] for parts 
of these expressions). The value of A^\m e jm^ 
obtained including these additional terms per- 
fectly agrees with that in Eq. (9) determined with 
the exact formulae. Indeed, their difference is of 
O(10~ 23 ), to be compared with the O(10~ 13 ) er- 
ror 8A^(m e /mf^) due to the present uncertainty 
of the ratio m e /m^. Therefore, it will be possible 
to compute (m e /m M ) with the simple expan- 
sion in Ref. [5] — thus avoiding the complexities of 
the exact expressions — even if the precision of the 
ratio m e /m^ will improve in the future by orders 
of magnitude. 

The contribution of the three-loop diagrams 
with both /i and t loop insertions in the photon 
propagator can be calculated numerically from 
the integral expressions of Ref. [17]. The value 
reported in Ref. [5] is 

A^^e/m^mjmr) = 1.909 45 (62) x 10~ 13 ,(11) 

a totally negligible O(10~ 21 ) contribution to 
a« ED . Adding up Eqs. (8), (9), (10) and (11) one 
obtains the three-loop QED coefficient [5] 

Ci 6) = 1.181234016827(19). (12) 

The relative contribution to a|? ED of the mass- 
dependent part of this three-loop coefficient is 
~0.1 ppb. This is smaller than the present 
~0.7ppb experimental uncertainty [1]. The error 
1.9xl0 -11 in Eq. (12) leads to a totally negligible 
O(10~ 19 ) uncertainty in a« ED . 

3.2. Determination of oc from a e 

As we already mentioned, recently a new mea- 
surement of the electron anomalous magnetic 
moment by Gabrielse and his collaborators 
achieved the extremely small relative uncertainty 
of 0.7 ppb [1], 

a EXP = 115 965 218 0.85 (76) x 10~ 12 . (13) 
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This uncertainty is nearly six times smaller than 
that of the last measurement of a e reported back 
in 1987, a EXP = 1159652188.3(4.2)xl(T 12 [30,10]. 
These two measurements differ by 1.7 standard 
deviations. 

The fine-structure constant a can be deter- 
mined equating the theoretical SM prediction of 
the electron g— 2 with its measured value 

a™(a)=a? p . (14) 

The SM prediction contains the QED contribu- 
tion a« ED (a) = ELiC^V/ 71 ") 1 (higher-order 
coefficients are assumed to be negligible), plus 
small weak and hadronic loop effects: al M (a) = 
a« ED (a) + a EW + < AD (the dependence on a of 
any contribution other than aJ? ED is negligible). 
The EW contribution is [10]: 

a EW = 0.0297(5) x 10~ 12 ; (15) 

this precise value includes the two-loop contribu- 
tions calculated in Refs. [31-33]. The hadronic 
term is [10,34,35]: 

< AD = 1.671(19) x 10- 12 . (16) 

The latest value for the four-loop QED coefficient 
is d S) = -1.7283(35) [36]. Following the argu- 
ment of [10], the educated guess ci 10) = 0.0(3.8) 
for the five-loop coefficient can be adopted. The 
errors 5Cf ] = 0.0035 and SCi W) = 3.8 lead to 
an uncertainty of 0.1 x 10~ 12 and 0.3 x 10~ 12 in 
a QED^ reS p ec tively. Solving Eq. (14) with the new 
measured value of Eq. (13), one obtains [4,5] 

aT 1 = 137.035 999 709(12) (30) (2) (90) 

= 137.035 999 709(96) [0.70 ppb]. (17) 

The first and second errors are due to the uncer- 
tainties of the four- and five-loop QED coefficient 
5Ce 8 ^ and 5Ci 10 \ respectively; the third one is 
caused by the tiny Sa" AD , and the last one (90 x 
10 -9 ) is from the experimental (5a EXP in Eq. (13). 
The uncertainties of the EW and two/three- 
loop QED contributions are totally negligible at 
present. 

The amazing precision of the determination in 
Eq. (17) represents the first significant improve- 
ment of this fundamental constant in a decade. 



At present, the best determinations of a inde- 
pendent of the electron g— 2 are 

aT 1 (Rb) = 137.035 998 78(91) [6.7 ppb], (18) 
aT 1 (Cs) = 137.036 000 0(11) [8.0 ppb]; (19) 

they are less precise by roughly a factor of ten. 
The value a _1 (Rb) was deduced from the mea- 
surement of the ratio h/M^b based on Bloch oscil- 
lations of Rb atoms in an optical lattice (h is the 
Planck constant and M-^h is the mass of the Rb 
atom) [37], while a _1 (Cs) was determined from 
the measurement of the ratio h/Mc s (Mc s is the 
mass of the Cs atom) via cesium recoil measure- 
ment techniques [38,39]. These two determina- 
tions of a also rely on the precisely known Ryd- 
berg constant and relative atomic masses of the 
electron, Rb and Cs atoms [10,40]. The values 
of a in Eqs. (18) and (19) are in good agreement 
with the result of Eq. (17), differing from the lat- 
ter by —1.0 and +0.3 standard deviations, respec- 
tively. This comparison provides a beautiful test 
of the validity of QED. It also probes for possible 
electron substructure [4]. 

4. MUON 

In this section we will summarize the present 
status of the SM prediction of the muon g— 2. See 
Rcfs. [41,42,21] for recent reviews. 

4.1. QED and Electroweak Contributions 

As we discussed in Sees. 2-3, the one-, two-, 
and three- loop QED contributions to a; (I = e, 
\i or t) are known analytically. The four-loop 
QED part of the muon 5— 2, which is about six 
times larger than the present experimental un- 
certainty of a^, has thus far been evaluated only 
numerically. This formidable task was first ac- 
complished by Kinoshita and his collaborators 
in the early 1980s [43]. The latest analyses ap- 
peared recently in [36,44]. The leading five- 
loop QED terms were recently evaluated in [45] 
and were presented at this workshop [46]; esti- 
mates obtained with the rcnormalization-group 
method agree with these results [47]. The to- 
tal QED prediction currently stands at a^ ED = 
116 584 718.09 (14) (08) x lO"' 11 [5]. The first error 
is determined by the uncertainties of the QED co- 
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efficients (dominated by the five- loop one), while 
the second is caused by the tiny uncertainty of the 
recent value of the fine-structure constant a [4,5]. 

The EW contribution to a M is suppressed by 
a factor (m^/M^) 2 with respect to the QED 
effects. The one-loop part was computed in 
f972 by several authors [48]: a™(l loop) = 

[l + |(1- 4sin 2 ^) 2 + 0(ml/Ml W:H ) , 

where G M = 1.16637(1) x 1CT 5 GcV~ 2 is 
the Fermi coupling constant and 9 W is the 
weak mixing angle. Closed analytic expres- 
sions for a™(l loop) taking exactly into ac- 
count the m^/Mj dependence (B = Z, W, 
Higgs, or other hypothetical bosons) can be 
found in Refs. [49]. Employing the on-shcll 
definition sin 2 6> w = 1 - M^/Mf [50], where 
M z = 91.1875(21) GcV and M w is the SM 
prediction of the W mass (which can be de- 
rived, for example, from the simple formu- 
lae of [51] leading to M w = 80.383 GeV for 
the Higgs mass M H = 150 GcV), one obtains 
a™ (1 loop) = 194.8 x 10" 11 . 

The two-loop EW contribution to a M is not neg- 
ligible because of large factors of ln(M z w /m/), 
where m/ is a fermion mass scale much smaller 
than M w [52]. It was computed in 1995 [31,32]. 
The proper treatment of the contribution of the 
light quarks was addressed in [53,54]. These re- 
finements significantly improved the reliability of 
the fcrmionic part (that containing closed fermion 
loops) of a® w (two- loop) leading, for M H = 
150 GeV, to a™ = 154(1) (2) x lO" 11 [54]. The 
first error is due to hadronic loop uncertainties, 
while the second one corresponds to an allowed 
range of M H e [114, 250] GeV, to the current top 
mass uncertainty, and to unknown three-loop ef- 
fects. The leading-logarithm three-loop contribu- 
tion to a™ is extremely small [54,55]. The re- 
sult of [56] for the two- loop bosonic part of a^ w , 
performed without the large M H approximation 
previously employed in Ref. [31], agrees with this 
previous evaluation in the large Higgs mass limit. 

4.2. The Hadronic Contribution 

The evaluation of the hadronic leading-order 
contribution a™°, due to the hadronic vacuum 
polarization correction to the one-loop diagram, 



involves long-distance QCD for which perturba- 
tion theory cannot be employed. However, using 
analyticity and unitarity, it was shown long ago 
that this term can be computed from hadronic 
e+e~ annihilation data via the dispersion inte- 
gral [57] 

1 1"°° 

<° = JZs dsK„(s)a(%), (20) 

where a^°\s) is the total cross section for e+e~ 
annihilation into any hadronic state, with extra- 
neous QED corrections subtracted off. The ker- 
nel function K^{s) decreases monotonically for 
increasing s. 

A prominent role among all e + e _ annihilation 
measurements is played by the precise data col- 
lected in 1994-95 by the CMD-2 detector at the 
VEPP-2M collider in Novosibirsk for the e + e~ — > 
7r + 7r~ cross section at values of yfs between 0.61 
and 0.96 GeV [58] (quoted systematic error 0.6%, 
dominated by the uncertainties in the radiative 
corrections). Recently [59] the CMD-2 Collab- 
oration released its 1996-98 measurements for 
the same cross section in the full energy range 
y/s e [0.37, 1.39] GeV. The part of these data for 
yfs G [0.61,0.96] GeV (quoted systematic error 
0.8%) agrees with their earlier result published 
in [58]. In 2005, also the SND Collaboration (at 
the VEPP-2M collider as well) released its analy- 
sis of the e+e~ — > 7r+7r~ process for yfs between 
0.39 and 0.98 GeV, with a systematic uncertainty 
of 1.3% (3.2%) for yfs~ larger (smaller) than 0.42 
GeV [60]. However, a recent reanalysis of these 
data [61] uncovered an error in the treatment of 
the radiative corrections, reducing the value of 
the measured cross section. The new SND result 
appears to be in good agreement with the corre- 
sponding one from CMD2. 

In 2004 the KLOE experiment at the DAFNE 
collider in Frascati presented a precise measure- 
ment of er(e + e~ — > n + n~) [62] via the initial- 
state radiation (ISR) method [63] at the <j> reso- 
nance. This cross section was extracted for y/s 
between 0.59 and 0.97 GeV with a systematic er- 
ror of 1.3% and a negligible statistical one. There 
are some discrepancies between the KLOE and 
CMD-2 results, even if their integrated contribu- 
tions to a" LO are similar. The data of KLOE [62] 
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and SND [61] particularly disagree above the p 
peak, where the latter are significantly higher. 
The study of the e + e~ — > 7r + 7r~ process via the 
ISR method is also in progress at BABAR [64] and 
Belle [65]. On the theoretical side, analyticity, 
unitarity and chiral symmetry provide strong con- 
straints for the pion form factor in the low-energy 
region [66]. Perhaps, also lattice QCD computa- 
tions of a^ LO , although not yet competitive with 
the precise results of the dispersive method, may 
eventually rival that precision [67]. 

Recent evaluations of the dispersive integral are 
in rather good agreement: 

a^ LO = 6909(39) e;cp (19) rad (7) 9cd xl0- 11 [68](21) 

a^ LO - 6894(42) e;cp (18) rad x KT 11 [69] (22) 

a™° = 6921 (56) x KT 11 [70] (23) 

a H ^° = 6944 (48) e;cp (10) rad x KT 11 [71]. (24) 

The three results in Eqs. (21)-(23) already in- 
clude the recent datasets from CMD-2 [59], 
SND [61] and BABAR [72] (not all of them, for 
the value in Eq. (23)). KLOE's data were not in- 
cluded in the analyses of Refs. [68] and [71]. The 
preliminary value in Eq. (21) updates the one of 
Ref. [73]. Significant progress is expected from 
the e + e~ collider VEPP-2000 under construction 
in Novosibirsk [74] and, possibly, from DAFNE-2 
at Frascati [75]. 

The authors of [76] pioneered the idea of using 
vector spectral functions derived from the study 
of hadronic r decays [77,78] to improve the eval- 
uation of the dispersive integral. However, the 
latest analysis with ALEPH, CLEO, and OPAL 
data yields a« LO = 7110 (50) e:cp (8) rQd (28) S[/(2) x 
10~ n [73], a value significantly higher than those 
obtained with e + e~ data. Isospin-breaking cor- 
rections were applied [79]. Although the pre- 
cise CMD-2 e + e~ — > n + n~ data are consistent 
with the corresponding r ones for energies be- 
low ~0.85 GeV, they are significantly lower for 
larger energies. KLOE's 7r + 7r~ spectral function 
confirms this discrepancy with the r data. Even 
if SND's 2005 results [60] were compatible with 
the r ones, their recent reanalysis [61], reducing 
the value of the measured cross section, seems to 
indicate that this is no longer the case. Interest- 
ing preliminary results from Belle of the study of 



the decay t~ — > ■k~tCv t were presented at this 
workshop [80,81]. This puzzling discrepancy be- 
tween the 7r + 7r~ spectral functions from e + e~ and 
isospin-breaking-corrected r data, also discussed 
at this workshop [82], could be caused by incon- 
sistencies in the e + e~ or t data, or in the isospin- 
breaking corrections which must be applied to the 
latter [83]. 

The hadronic higher-order (a 3 ) contribution 
a™° can be divided into two parts: a" HO = 
a™°(vp) + a« HO (lbl). The first one is the 0(a 3 ) 
contribution of diagrams containing hadronic vac- 
uum polarization insertions [35]. Its latest value 
is a£ HO (vp) = -97.9 (0.9) e;cp (0.3) rad x KT 11 [84, 
69]; it changes by <~ — 3 x 10 -11 if hadronic r- 
decay data are used instead [42]. The second 
term, also of 0(a 3 ), is the hadronic light-by- 
light contribution. As it cannot be directly de- 
termined via a dispersion relation approach us- 
ing data (unlike the hadronic vacuum polariza- 
tion contribution) , its evaluation relics on specific 
models of low-energy hadronic interactions with 
electromagnetic currents. Three major compo- 
nents of o™°(lbl) can be identified: charged-pion 
loops, quark loops, and pseudoscalar (ir°, 77, and 
rf) pole diagrams [85]. The latter ones dominate 
the final result and require information on the 
electromagnetic form factors of the pseudoscalars. 
In 2001 the authors of [86] uncovered a sign error 
in earlier evaluations of the dominating pion-pole 
part. Their estimate of a™°(lbl), based also on 
the previous results for the quark and charged- 
pions loop parts, is a™°(lbl) = 80 (40) x 10" 11 . A 
higher value was obtained in 2003 including short- 
distance QCD constraints: a™°(lbl) = 136 (25) x 
lO" 11 [87]. The upper bounds found in Refs. [88] 
are consistent with this higher value. Further in- 
dependent calculations would provide an impor- 
tant check of these results for a™°(lbl) (see [89] 
for a recent critical comparison of these evalua- 
tions), a contribution whose uncertainty may be- 
come the ultimate limitation of the SM prediction 
of the muon g— 2. 

4.3. Standard Model vs. Measurement 

The first column of Table 1 shows a^ M = 
a« ED + a™ + a* LO + a™°. The values employed 
for a™° are indicated by the reference on the 
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Table 1 



a^: Standard Model vs. measurement. 





x 10 11 


A x 10 11 


(7 


[68] 


116 591763 (60) 


317 (87) 


3.7 (3.2) 


[69] 


116 591748 (61) 


332 (88) 


3.8 (3.4) 


[70] 


116 591775 (69) 


305 (93) 


3.3 (2.8) 


[71] 


116 591798 (63) 


282 (89) 


3.2 (2.7) 


[73] 


116591961 (70) 


119 (95) 


1.3 (0.7) 



left; the value of Ref. [73] was obtained using 
data from hadronic t decays. All a^ M values were 
derived with o™°(lbl) = 80 (40) x lO" 11 [85,86]. 
Errors were added in quadrature. The present 
world average experimental value is a EXP = 
116 592 080(63) x lO" 11 [0.5 ppm] [3]. Note that 
the theoretical error on is now roughly the 
same as the experimental one. The differences 
A = a EXP — a^, M arc listed in the second column 
of Table 1, while the numbers of "standard devi- 
ations" (a) appear in the third one. Lower dis- 
crepancies, shown in angle brackets, are obtained 
if a™°(lbl) = 136 (25) x 10~ n [87] is used instead 
of 80(40) x 10~ n [85,86]. 

5. TAU 

This section summarizes the results reported at 
this workshop [90] for the theoretical prediction 
of the t lepton g— 2. We also refer the reader to 
Ref. [91] for a very recent update and review. 

5.1. QED and Electroweak Contributions 

The sum of the one-, two- and three-loop QED 
contribution to a T , obtained using the mass ratios 
given in Sec. 2 is [5]: 

a^ ED = 117324(2) x 10~ 8 . (25) 

The error da® ED is the uncertainty 8C^ (a/7r) 4 ~ 
7T 2 \n 2 (m T /m e )(a/n) 4 <~ 2 x 10~ 8 assigned in 
[5] to a^ ED for uncalculated four-loop contribu- 
tions. The errors due to the uncertainties of the 
0{a 2 ) (5 x 10- 10 ) and 0(a 3 ) terms (3 x 10" 11 ), 
as well as that induced by the uncertainty of a 
(8 x 10~ 13 ) are negligible. The result in Eq. (25) 
supersedes the earlier values a® ED = 117 319 (1) x 



10~ 8 [92] and a? ED = 117 327.1 (1.2) x 10~ 8 [93] 
(see Refs. [5,91] for details). 

The EW correction to the anomalous magnetic 
moment of the r lepton is of the same order of 
magnitude as the three-loop QED one. The one- 
loop EW term is a EW (l loop) = 55.1 x 10~ 8 [48]. 
The estimate of the total EW contribution of 
Ref. [92], a EW = 55.60(2) x 10~ 8 , obtained from 
the one-loop formula, is similar to the one-loop 
value reported above. However, it doesn't contain 
the two-loop contribution, which is not negligible. 

As we already discussed for the muon g— 2, the 
two-loop EW contributions af w (2 loop) (l = e, \i 
or r) were computed in 1995 [31-33], leading to 
a significant reduction of the one-loop prediction 
because of large factors of ln(M ZiW -/m/), where 
nif is a fcrmion mass scale much smaller than 
M w [52]. The numerical value of the bosonic 
part of a EW (2 loop), for M H = 150 GeV, is 
a EW (2 loop bos) = -3.06 x 10~ 8 [31,33,90,91]. 
Very recently, the analysis of the fermionic part 
of Ref. [32,33] was slightly refined in Ref. [91] 
leading, for M H = 150 GeV, to a EW (2 loop fcr) = 
—4.68 x 10~ 8 . The sum of the fermionic 
and bosonic two-loop EW contributions gives 
a EW (2 loop) = -7.74 x 10~ 8 [90,91], a 14% re- 
duction of the one- loop result. As discussed in 
Sec. 4, the leading-logarithm three-loop EW con- 
tributions to the muon g— 2 were determined to be 
extremely small via renormalization-group anal- 
yses [54,55]; in Ref. [90,91] an additional un- 
certainty of 0[a EW (2 loop)(a/7r)ln(Mj/m2)] ~ 
O(10~ 9 ) was assigned to a EW to account for these 
neglected three- loop effects. Adding a EW (2 loop) 
to the one-loop value presented above, one gets 
the total EW correction (M ff = 150 GeV) [90,91]: 

fl Ew = 47 4(5) x 10 -8 ^6) 

The uncertainty allows M H to range from 114 
GeV up to ^300 GeV, and reflects the estimated 
errors induced by hadronic loop effects, neglected 
two-loop bosonic terms and the missing three- 
loop contribution. It also includes the tiny errors 
due to the uncertainties in M top and m T . The 
value in Eq. (26) is in agreement with the predic- 
tion a EW = 47(1) x 10~ 8 [33,93], with a reduced 
uncertainty. As we mentioned earlier, the EW es- 
timate of Ref. [92], a EW = 55.60(2) x 10~ 8 , mainly 



8 



M. Passera 



differs from Eq. (26) in that it doesn't include the 
two-loop corrections. 

5.2. The Hadronic Contribution 

Similarly to the case of the muon, the leading- 
order hadronic contribution to the r lepton g—2 is 
given by the dispersion integral in Eq. (20) (with 
the kernel K fl (s) replaced by K T (s)), in which 
the role of the low energies is very important, al- 
though not as strongly as in a" LO . The history of 
the a" LO calculations [91-97] based mainly on ex- 
perimental e + e _ data is shown in Table 2. Purely 
theoretical estimates somewhat undervalue the 
hadronic contribution and have rather large un- 
certainties [98-101]. 



Table 2 




Calculations of a" L °. 




Author 


a« LO x 10 8 


Narison [94] 


370 ± 40 


Barish & Stroynowski [95] 


-350 


Samuel et al. [92] 


360 ± 32 


Eidelman & Jegerlehner [96,97] 


338.4 ±2.0 ±9.1 


Narison [93] 


353.6 ± 4.0 


Eidelman & Passera [90,91] 


337.5 ± 3.7 



The calculation of the leading-order contribu- 
tion was very recently updated using the whole 
bulk of experimental data below 12 GeV, which 
include, as discussed in Sec. 4, old data compiled 
in Refs. [96,73], and recent results from CMD- 
2 [59], SND [61], KLOE [62] and BABAR [72]. The 
improvement is particularly strong in the chan- 
nel e + e~ — > 7r + 7r~. The result of this analysis 
is [90,91]: 

< LO = 337.5(3.7) x 10~ 8 . (27) 

The overall uncertainty is 2.5 times smaller than 
that of the previous data-based prediction [96,97]. 

Like a™°, the hadronic higher-order (a 3 ) con- 
tribution a™° can be divided into two parts: 

fl HHO = a HH 0(vp ) + a HHO( lb l) ^ ^ 

the 0(a 3 ) contribution of diagrams containing 
hadronic self-energy insertions in the photon 
propagators, is [35]: 

a™°(vp) = 7.6(2) xl0~ 8 . (28) 



Note that naively rescaling the muon result by 
the factor m 2 /m 2 l leads to the totally incorrect 
estimate a" HO (vp) = (-101 x 10~ n ) x m 2 T /m 2 ^ = 
-29 x 10~ 8 (the a™°(vp) value is from Ref. [35]); 
even the sign is wrong! Until recently, very 
few estimates of the light-by-light contribution 
a™°(lbl) existed in the literature [92,93,35], and 
all of them were obtained simply rescaling the 
muon results a™° (lbl) by a factor m 2 / m 2 ^ . These 
very naive estimates fall short of what is needed, 
as this scaling is not justified. For these reasons, 
a parton- level estimate of a™°(lbl) was recently 
performed in Ref. [91], obtaining 

a™°(lbl) = 5(3) x 10~ 8 . (29) 

This value is much lower than those obtained by 
simple rescaling. We refer the reader to Ref. [91] 
for a detailed discussion. 

5.3. Standard Model prediction for a T 

We can now add up all the contributions dis- 
cussed in the previous sections to derive the SM 
prediction for a T [90,91]: 

fl SM = a QED + fl EW + fl HLO + ffl HHO 

= 117721 (5) x 10~ 8 . (30) 

Errors were added in quadrature. Comparing 
the most stringent experimental limit mentioned 
above, -0.052 < o? xp < 0.013 at 95% confidence 
level [7], with Eq. (30), it is clear that the sensitiv- 
ity of the best existing measurements is still more 
than an order of magnitude worse than needed. 
For other limits on a T see Ref. [13,102]. 

6. CONCLUSIONS AND OUTLOOK 

The study of the anomalous magnetic moment 
of the electron and the muon provides a beauti- 
ful example of interplay between theory and ex- 
periment. It started more than 50 years ago: 
Schwingcr's 1948 calculation [9] of the leading 
contribution to a e was one of the very first results 
of QED, and its agreement with the experimental 
value provided one of the early confirmations of 
this theory. In July 2006 Gabrielse and his collab- 
orators published their measurement of a e with 
an unprecedented accuracy of 0.7 ppb [1], break- 
ing the 1987 world record held by Dehmelt and his 
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collaborators [30] . On the theoretical side, heroic 
higher-order QED calculations were performed 
over the years by Kinoshita, Remiddi, and their 
collaborators. A result of this extraordinary syn- 
ergy is today's determination of the fine-structure 
constant via o e , a" 1 = 137.035 999 709(96) [4,5], 
whose uncertainty is roughly ten times smaller 
than that from any other method. A further re- 
duction of the experimental error <5<2g XP , which 
dominates the uncertainty of a, appears to be at 
hand [1,4]. 

The SM prediction of the muon g—2 deviates 
from the present experimental value by more 
than 3cr, if data from e + e~ collisions are em- 
ployed to evaluate the leading-order hadronic 
term (the recent datasets from the CMD-2 [59], 
SND [61] and BABAR [72] are already included). 
This deviation cannot be ignored. Determina- 
tions based on r-decay data deviate only by ap- 
proximately 1 a. The puzzling discrepancy be- 
tween the tt + tt~ spectral functions from e + e~ 
and isospin-breaking-corrected r data could be 
caused by inconsistencies in the e + e~ or r data, 
or in the isospin-breaking corrections applied to 
the latter. Indeed, some disagreements occur be- 
tween e+e - data sets, requiring further detailed 
investigations. On the other hand, the connection 
of t data with the leading hadronic contribution 
to a M is less direct, and one wonders whether all 
possible isospin-breaking effects have been prop- 
erly taken into account. 

The impressive results of the E821 experiment 
are still limited by statistical errors. A new ex- 
periment, E969, has been approved (but not yet 
funded) at Brookhaven in 2004 [103] . Its goal is to 
reduce the present experimental uncertainty by a 
factor of 2.5 to about 0.2 ppm. A letter of intent 
for an even more precise muon g—2 experiment 
was submitted to J-PARC with the proposal to 
reach a precision below 0.1 ppm [104]. While the 
QED and EW contributions appear to be ready to 
rival these precisions, much effort will be needed 
to reduce the hadronic uncertainty by a factor of 
two. This effort is challenging but possible, and 
certainly well motivated by the excellent oppor- 
tunity the muon g—2 is providing us to unveil (or 
constrain) NP effects. 

The g—2 of the r lepton is even more sensitive 



than the muon one to EW and NP loop effects 
that give contributions ^mf. However, unfor- 
tunately, the very short lifetime of the r lepton 
makes it very difficult to determine a T by measur- 
ing its spin precession in the magnetic field, like 
in the muon g—2 experiment [2,3]. Instead, ex- 
periments focus on high-precision measurements 
of the r lepton pair production in various high- 
energy processes and comparison of the measured 
cross sections with the QED predictions [7,13], 
but their sensitivity is still more than an order 
of magnitude worse than that required to deter- 
mine a T . 

Nonetheless, there are many interesting sugges- 
tions to measure a T , e.g., from the radiation am- 
plitude zero in radiative r decays [105] or from 
other observables; these methods could possibly 
exploit the very large r lepton samples collected 
at B factories. Reference [106] suggests a simi- 
lar method to study a T using radiative W decays 
and potentially very high data samples at LHC. 
Yet another method would use the channeling in 
a bent crystal similarly to the measurement of 
magnetic moments of short-living baryons [107]. 
In the case of the r lepton, it was suggested to 
use the decay B + — > t + v t , which would pro- 
duce polarized r leptons [92] and was recently ob- 
served [108]. We believe that a detailed feasibility 
study of such experiments, as well as further at- 
tempts to improve the accuracy of the theoretical 
prediction for a T , are quite timely. 
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